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ABSTRACT 



Context. Nearby edge-on galaxies showing a synchrotron halo are nearly ideal objects for studying the transport of cosmic rays (CRs) 
in galaxies. Among them, the nearby starburst galaxy NGC 253 hosts a galactic wind indicated by various ISM phases in its halo. 
Aims. The diffusive and convective CR transport from the disk into the halo is investigated using the local CR bulk speed. The 
connection between the CR transport and the galactic wind is outlined. 

Methods. We observed NGC 253 with the VLA at /16.2 cm in a mosaic with 15 pointings. The missing zero-spacing flux density of 
the VLA mosaic was filled in using observations with the 100-m Effelsberg telescope. We also obtained a new /13.6 cm map from 
Effelsberg observations and reproduced VLA maps at A20 cm and /190 cm. The high dynamic range needed due to the strong nuclear 
point-like source was addressed with a special data calibration scheme for both the single-dish and the interferometric observations. 
Results. We find a thin and a thick radio disk with exponential scaleheights of 0.3 kpc and 1 .7 kpc at /16.2 cm. The equipartition total 
magnetic field strength between 7 jiG and 1 8 //G in the disk is remarkably high. We use the spectral aging of the cosmic ray electrons 
(CREs) seen in the vertical profiles of the spectral index to determine a lower limit for the global CR bulk speed as 170 ± 70 km s -1 . 
The linear correlation between the scaleheights and the CRE lifetimes, as evident from the dumbbell shaped halo, requires a vertical 
CR transport with a bulk speed of 300 ± 30kms~' in the northeastern halo, similar to the escape velocity of 280kms~'. This shows 
the presence of a "disk wind" in NGC 253. In the southwestern halo, the transport is mainly diffusive with a diffusion coefficient of 
2.0 ± 0.2 x 10 29 cm 2 s _1 . 

Conclusions. Measuring the radio synchrotron scaleheight and estimating the CRE lifetime allow us to determine the bulk speed 
of the CR transport into the halo. The transport is convective and more efficient in the northeastern halo, while it is diffusive in the 
southwestern halo. The luminous material is transported by the disk wind, which can explain the different amounts of extra-planar 
HI, Ha, and soft X-ray emission in the two halo parts. Future low-frequency radio observations will provide the data to analyze the 
vertical velocity profile of galactic winds. 
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1. Introduction 

NGC 253 i s a prototypical starburs t galaxy at a distance of 
3.94 Mpc dKarachentsev et all I2003I) . It is a late-type spiral 
galax y with a high inclination angle of 78.5° ± 0.5° dPencd 
1980). This allows observations of extra-planar components of 
the interstellar medium (ISM) at different wavelengths: ROSAT 
Position Sensitive Proportional Counter (PSPC) observations 
have revealed a very extended halo of soft diffuse X-ray emis- 
sion extending into the halo to a distance of 8 kpc from the 
disk. The prominent lobes of soft X-ray emissio n show a dis- 
tinct "X"-shaped pattern centered on the nucleus (Pietsch et all 



|2000|) . Moreover, Ha observations bv lHoopes et al. ( 19961) show 
plumes emerging from the disk into the halo, which spatially co- 
incide with the lobes o f the soft X-ray emission. Observations of 
iBoomsma et al.l d2005l) in H I show large plumes of extra-planar 
H I emission enclosing the soft X-ray emission. 

Given the different phases of the ISM in the halo of 
NGC 253, the question arises of how the observed halo struc- 
ture was formed. The distribution of the different ISM phases in 
shells has led to the suggestion that the halo can be considered as 
one huge bubble that is fueled by an active star-forming region 
around the nucleus of the galaxy ([Strickland et all 12002b . This 



scenario is supported by the finding of a small plume of hot X- 
ray gas in CHANDRA m aps associated with a conical outflow 
from the nuclear starburst dStrickland et al.l200 0). Spectroscopic 
measurement s of Ha-emitting gas in the southern nuclear out- 
flow cone bv lSchulz & Wegnerl dl992^ give an outflow vel ocity 
of 390 km s _1 . Numerical models by ISuchkov et al.l d 19941) fol- 
low the evolution of a multi-phase halo powered by a compact 
source located in the center of the disk. A close resemblance of 
the "X"-shaped pattern of the soft X-ray emission was indeed 
found, which could be interpreted as the termination shock of 
the hot gas where the gas cools effectively. 

A possible explanation for the amount of X-ray and Ha- 
emitting gas in th e halo is the dissipation o f kinetic energy 
via shock heating dLehnert & Hec kman 1996) if the ther mal- 
ization fraction is high enough dStrickland et al.l 120021) . A 
precondition for this model is the existence of cool gas in the 
halo that condenses to clouds from whi ch the shock ca n form . 
ISOPHOT observations discussed by iRadovich etafj (1200 ll) 
show "X"-shaped far-infrared emission overlapping with the 
Ha and X-ray emission, indicating the existence of such a 
cool gas. However, the origin of the gas is not clear. It could 
be a preexisting population of cool gas in the halo, possibly 
accreted from a merger with a satellite. Another possibility 
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is that the gas of the disk can be transported into the halo 
along the walls of the cavity as Kelvin-Helmho ltz instabilities 
rip off clumps of cool gas in the disk (see e.g. ISuchkov et alj 
fl99l iHeckman et all 12000). Summarizing, it is important to 
understand how the gas is transported from the disk into the halo. 

Cosmic rays (CRs) and magnetic fields are important for the 
evolution of the ISM, because they may poss ess an energy den - 
sity comparable to that of the thermal gas dBeck et alj 1 1996h . 
They are even more important for the study of galactic winds, be- 
cause the thermal gas itself can drive a galactic wind only when 
the heating processes are sufficiently strong and the gravitational 
potential is comparably weak. A good example of this situa- 
tion is the edge-on galaxy M 82, which possesses a huge halo 
of hot X-ray, Ha, and radio continuum emission. However, the 
hot gas cools radiatively preventing the escape from the gravita- 
tional potential and falls back onto the disk in condensed clouds . 
This galactic fountain scenario proposed by iField et ail d 19691) 
might explain the infall of the so-called High Velocity Clouds 
that could be interpreted as cooled remnants of the hot gas. 

If we incorporate the CR pressure the galactic wi nd can be 
sustained against the radiative losses dlpavichlll975l) . because 
the CRs suffer only little from radiative losses (mainly via syn- 
chrotron losses of the electrons). As the protons carry most of the 
energy this does not influence the CR energy density. Thus the 
pressure support from below can drive a g alactic wind against 
the gr a vity o f the underlying galactic disk. iBreitschwerdt et all 
d 1 99 It 1 19931) started a series of papers treating the CR trans- 
port along "flux density tubes" in a one-dimensional model. This 
model was extended to a full three-dimensional t reatment (with 
rotational symmetry) by Zir akashvili et all ( 1 19961) . including the 
influence of galactic rotation. It was found that CRs can drive 
a galactic wind under a broad variety of conditions in the disk 
and that the loss of m ass and momentum over the lifetime of 
a galaxy is s i gnifica nt (iPtuskin et al.lll997l) . A recent study by 
lEverett et alj (|2008) showed that the extra-planar X-ray emis- 
sion in the Milky Way can be explained by a CR driven galactic 
wind. 

The theory of a CR driven galactic wind can be constrained 
with radio continuum studies: as the relativistic CRs are the most 
pervasive component of the ISM, they are expected to be associ- 
ated with the other phases in the halo of NGC 253. Indeed, a ra- 
dio halo around NGC 253 was detected bv lBeck et ail d 19791) us- 
ing the 100-m EfFelsberg telescope at A3. 6 cm. However, no cor- 
rection for the sidelobes of the strong nuclear point-like source 
was applied and thus it was not possible to disentangle the in- 
trinsic halo emission from the contribution by the sidelobes. The 
st udy was extended la ter with VLA interferometric observations 
bv lCarilli et al.l d 1992b . which show an extended radio halo over 
the length of the optical disk very different in shape from the halo 
seen in X-rays. The finding of a large "radio spur" southeast of 
the nucleus already led to speculations of an outflow connected 
with the spur. 

Extending the observations to other wavelengths and includ- 
ing polarization, a predict ion of the magnet ic field structure in 
the disk was presented bv lBeck et al.1 d!994l) . The sensitivity of 
the observations, however, was not sufficient to get information 
for the extra-planar polarized emission so that the magnetic field 
structure in the halo could not be determined. With the new VLA 
mosaic at A6.2 cm we are able to increase the sensitivity at short 
wavelengths, allowing us to study in detail the extra-planar emis- 
sion and possible relations to the outflow of hot gas. The ad- 
vantage is that the high synchrotron losses at this short wave- 
length impose tighter constraints on the CR transport from the 



disk into the halo. Moreover, for the first time sensitive polar- 
ization measurements were carried out to study the large-scale 
magnetic field structure in the halo of NGC 253, which will be 
addressed in a subsequent paper (thereafter Paper II). 

This paper is organized as follows: in Sect. [2] we introduce 
our observations and present the data reduction. We present our 
radio continuum maps along with a discussion of the morphol- 
ogy in Sect. [3] The distribution of the radio spectral index is 
presented in Sect. [4] The finding of an extended radio halo with 
a thick radio disk leads to the determination of the vertical CR 
bulk speed (Sect. |5). In Sect. [6] we discuss the impact of our 
findings on the feeding of the radio halo with CRs from the star- 
forming disk. Finally, we summarize our results and present our 
conclusions in Sect. [7] 

2. Observations and data reduction 

2.1. The Effelsberg observations 

For the single-dish observations at /13.6 cm and /16.2 cm of 
NGC 253 we used the Effelsberg 100-m telescopeQ . The A3.6 cm 
observations were obtained between February 2003 and July 
2004 with a single-horn receiver mounted in the secondary focus 
cabin (SFK). We used the broad band receiver channel centered 
at 8.35GHz which has a bandwidth of 1.1 GHz. The /16.2cm 
observations were carried out in 1997 with a dual-horn receiver 
mounted in the SFK using a bandwidth of 500 MHz centered at 
4.85 GHz. Both receivers are capable to record the total intensity 
and the Stokes channels Q and U simultaneously. As primary 
(fl ux) calibrator we u sed 3C286 adopting the flux density scale 
of iBaars et al 1 (119771) . 

We observed 26 maps at /13.6 cm with a map size of 30' x 
20'. Two perpendicular scanning directions were aligned with 
the major / minor axis of NGC 253 (p.a. = 52°), in order to 
reduce the time for scanning. For the A 6.2cm observations we 
m ade use of the "softwar e beam switching" technique described 
in iMorsi & Reichl (Qj)86). Since the observations with the dual- 
horn are made in the azimuth- elevation coordinate s ystem, we 
used a restoration method (see lEmerson et al. I [19791) to obtain 
eleven maps with a size of 30' x 30'. 

The data r eduction was p erformed using NOD2 data reduc- 
tion package tH aslamlll974l) . Each map was checked and cor- 
rected for pointing errors using Gaussian fits on the nuclear 
point-like source in NGC 253. The final map was produced using 
the P LAIT algorithm, which w e used to suppress scanning ef- 
fects dEmerson & GraveJI 19881) . The noise level of the final map 
at /13.6 cm is 500 /iJy /beam in total power with a resolution of 
84" 0. At /16.2cm a noise level of 1 mJy/beam with a resolution 
of 144" was achieved. 

The sidelobes of the Effelsberg beam pattern peak at a level 
of a few % which limits the dynamic range to 100 at most. As 
the nuclear point-like source has a flux density of 1 Jy, we would 
be limited to a noise level of 10 mJy/beam. In order to overcome 
these limitations we had to remove the sidelobes of the nuclear 
point-like source. This was accomplishe d by applying a p roce- 
dure similar to the CLEAN algorithm bv lHogboml dl974l) . This 
procedure was revised and extended for applica tion on single- 
dish observations by K lein and Rottmann (see iKlein & Mackl 
1 19951 : lRottmannlll996l) . The required beam patterns were ob- 
tained from deep observations of 3C84 and 0234-23 at /13.6 cm 

1 The Effelsberg 100-m telescope is operated by the Max-Planck 
Institut fur Radioastronomie (MPIfR). 

2 All angular resolutions in this paper are referred to as the Half 
Power Beam Width (HPBW). 
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Fig. 1. The VLA mosaic of NGC 253 at /16.2 cm consisting of 
15 pointings. The circles indicate the size of the VLA primary 
beam which is 8' at A6.2 cm. The background optical image is 
from the DSS. 



and /16.2 cm, respectively. Each single map of NGC 253 was 
cleaned with the mean parallactic angle of the map, respectively. 
This is sufficient, since the variation of the parallactic angle is 
less than 5° during the observation of a single map. 

We were able to obtain essentially noise limited maps with 
a dynamic range of at least 1000. The cleaned structures (the 
difference between the uncleaned and the cleaned map) exhibit 
prominent sidelobes extending far from the disk underlining the 
importance of a careful cleani ng procedure. D etails of the single- 
dish cleaning can be found in lHeeserJ (120081) . 

2.2. The VLA observations 

With the VLAEI in D-configuration we observed NGC 253 in July 
2004 at /16.2 cm with a total integration time of 40 hours. The 
size of the primary beam of the VLA antennas at A6.2 cm is 8', 
requiring a mosaic of 15 pointings. They were arranged in three 
rows parallel to the major axis of the galaxy, including much of 
the extra-planar field of interest. An illustration of the mosaic 
setup is shown in Fig. Q] as an overlay onto an optical image 
from the DSS^ The data were calibrated with 3C48 and 3C138 
as primary (flux) calibrators using the standard procedures in the 
Astronomica l Image Processing System^ with the flux density 
scale given in lBaars et al.ldl977l) . The performed self-calibration 
allowed us to suppress the sidelobes of the nuclear point-like 
source in order to gain a high dynamic range. 



3 The VLA (Very Large Array) is operated by the NRAO (National 
Radio Astronomy Observatory). 

4 The compressed files of the "Palomar Observatory - Space 
Telescope Science Institute Digital Sky Survey" of the northern sky, 
based on scans of the Second Palomar Sky Survey are ©1993-1995 
by the California Institute of Technology and are distributed herein by 
agreement. 

3 The Astronomical Image Processing System (AIPS) is distributed 
by the National Radio Astronomy Observatory (NRAO) as free soft- 
ware. 
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Fig. 2. Total power radio continuum at /16.2 cm obtained from 
the VLA observations with 30" resolution. Contours are at 3, 
6, 12, 24, 48, 96, 192, 384, 768, 1536, 3077, 6144, 12288, and 
24576 x 30 yuJy/bea m. The grey scale shows Ha emission from 
iHoopes et all dl996l) . 
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Fig. 3. Total power radio continuum at /16.2 cm obtained from 
the Effelsberg observations with 144" resolution. Contours are 
at 3, 6, 12, 24, 48, 96, 192, 384, 768, and 1536 x 1 mJy/beam. 
The grey scale shows Ha emission as in Fig. [2] 



After inverting the (u,v)-data and the successive cleaning we 
convolved the maps of each pointing with a Gaussian in order 
to obtain an identical clean beam. Using LTESS (part of AIPS) 
the maps were finally combined with a linear superposition and a 
correction for the VLA primary beam attenuation using informa- 
tion from eac h pointing out to the 7 % level of the primary beam 
(lBraunlll988l) . Our final map has a noise level of 30yuJy/beam 
in total power which corresponds to a dynamic range of 40000 
with respect to the nucleus. 
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Fig. 4. Total power radio continuum at /(90 cm obtained from 
VLA observations with 70" resolution. Contours are at 3, 6, 12, 
24, 48, 96, 192, 384, 7 68, 1536, and 3077 x 8mJy/beam. The 
map is reproduced from lCarilli et al.l (119921) . 
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Fig. 5. Total power radio continuum at /120 cm obtained from 
VLA observations with 30" resolution. Contours are at 3, 6, 12, 
24, 48, 96, 192, 384, 76 8 1536, and 3077x 0.35 mJy/beam. The 
map is reproduced from lCarilli et al.l d!992l) . 
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Fig. 6. Total power radio continuum at /16.2 cm obtained from 
combined VLA and Effelsberg observations with 30" resolution. 
Contours are at 3, 6, 12, 24, 48, 96, 192, 384, 768, 1536, 3077, 
6144, 12288, and 24576 x 30/Jy/beam. 



An important shortcoming of interferometric observations 
is the limited sensitivity to extended emission on large angular 
scales. The VLA in D-configuration at /16.2 cm is only sensi- 
tive to emission with a scale smaller than 5'. Since NGC 253 is 
much larger, we observe with the VLA only a fraction of the 
extended emission detected with a single-dish telescope. The in- 
fluence of the missing zero-spacing flux density can be easily 
seen by comparing the Effelsberg map at A6.2 cm in Fig. [3] with 
the VLA map shown in Fig. [2] The missing total flux density is 
about 10% of the total flux density measured with Effelsberg; 
if the flux density of the extended emission (without the nuclear 



Fig. 7. Total power radio continuum at /13.6 cm obtained from 
Effelsberg observations with 84" resolution. Contours are at 3, 
6, 12, 24, 48, 96, 192, 384, 768, and 1536 x 0.5 mJy/beam. 



point-like source) is considered, the fraction of the missing zero- 
spacing flux density increases to 20 %. 

We used the Effelsberg map in order to fill in the missing 
zero-spacing flux in the VLA map. This was done using IMERG 
(part of AIPS) which performs a Fourier combination of the two 
maps. We stress that in particular for mosaics, which cover the 
entire region of interest, a single-dish map is necessary in order 
to detect all extended emission. 

The /190 cm and A 2Q cm maps are from VLA observations 
bv lCarilli et al.l d 1992b . While one pointing with a primary beam 
size of 138' is sufficient to cover NGC 253 at /190 cm, three 
pointings with a beam size of 25' were combined to a mosaic 
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at A20 cm. The largest angular scale structure which can be im- 
aged reasonable well with the VLA is 70' at /190 cm and 15' 
at A2Q cm, respectively. The needed correction for the missing 
zero-spacing flux of the mosaic at A2Q cm was not performed be- 
cause no single-dish map is available at this wavelength. Hence, 
the A20 cm map is not well suited to investigate the structure of 
the extended emission in the halo. 



3. Morphology of the total power emission 

We present the total power maps as contour plots with the first 
contour indicating emission at three times the rms noise level 
(3 rms x 2 m_1 ; m = 1, 2, 3, . . .). The optical background image 
was obtained from the DSS. A summary of the maps presented 
in this Paper can be found in Table Q] 



Table 1. The maps of NGC 253 included in this paper. 



A 

[cm] 


resolution 


rms 

[mJy/beam] 


dynamic 
range 


references 




90 


70" 


8 


500 


Carilli et al. 


(1992) 


20 


30" 


0.35 


5700 


Carilli et al. 


(1992) 


6.2 


30" 


0.03 


42000 


this paper 




3.6 


84" 


0.5 


2000 


this paper 





The continuum emission not only extends over the optical 
disk along the major axis, but we also find a notable exten- 
sion of emission along the minor axis in all maps. We refer to 
this as extra-planar emission. A detailed analysis and proof that 
the emission does not belong to the disk is given in Sect. 15.11 
The extent of extra-planar emission is strongly dependent on the 
observation wavelength, which can be seen by comparing the 
maps. The /190 cm map shown in Fig. |4]has very pronounced 
extra-planar emission, while the disk does not show up promi- 
nently. The A2Q cm and A6.2 cm maps presented in Figs. [5] and 
[6] have less extra-planar emission and at A3. 6 cm there is almost 
no extra-planar emission (Fig. [7}. 

The morphology of the total power emission also depends 
on the resolution and the sensitivity of the maps. In the /120 cm 
and A6.2 cm maps with a resolution of 30" the spiral arms are 
distincted features in the disk, one extending to the northeast and 
another to the southwest. The inner part of the disk is visible as a 
plateau of high emission coinciding with the brightest part in the 
optical image. In all maps, independent of the resolution, there is 
an extension southeast of the nucleus with an unresolved source 
at its southern tip. This is an example of a so-called "radio spur", 
where the magnetic field has a strong vertical component, as will 
be shown in Paper II. 

The total extent of emission perpendicular to the disk is not 
a good measure for the amount of extra-planar emission since 
it depends on the sensitivity of the observations. This is visi- 
ble in the A2Q cm map which shows less extra-planar emission 
than the A6.2 cm map because its noise level is 10 times higher. 
Thus, with a low sensitivity we observe only the brightest tip 
of the radio halo. Moreover, the extent appears to be larger for 
low resolutions. This is the motivation to use the scaleheight h 
which prescribes the emission profiles perpendicular to the disk 
as an exponential function oc exp(-z//i) of the distance z from 
the galactic midplane. In Sect. l5.ll we discuss the scaleheights in 
order to characterize the radio halo. 



We find that the extent of the extra-planar emission varies 
with the offset along the major axis. The contour lines in all 
maps show a dumbbell shape which means that the extension 
of the extra-planar emission has a local minimum at the center 
of the disk. While this trend is visible in all maps it is clearly 
seen at the shorter wavelengths. Especially the A3. 6 cm map has 
a prominent neck above and below the nucleus, where the con- 
tour lines indicate a steep gradient along the minor axis. This 
points to a smaller scaleheight of the extra-planar emission. If 
the scaleheight were constant along the major axis, the largest 
vertical extension of the emission would be expected in the cen- 
ter. We investigate this issue and give an explanation in Sect. [5] 

4. Spectral index distribution 

4.1. Unresolved sources 

Our continuum maps, in particular those with the highest resolu- 
tion, show several unresolved (point-like) sources. We identified 
the sources in order to subtract them from the maps prior to our 
analyis of the extended emission. We obtained peak flux densi- 
ties and positions by fitting Gaussians to them, where we cor- 
rected the peak flux density for the underlying extended emis- 
sion. This is particularly important for the nuclear point-like 
source, because it is embedded in strong extended emission. The 
list of point-like sources presented in Table[2]includes positions, 
flux densities, and spectral indices. The spectral a is defined by 
the relation S oc v~ a , where S is the total power radio continuum 
flux density and v is the observation frequency. The strongest 
two sources, the nucleus (no. 3) and the source in the radio 
spur (no. 6), are detected at all four wavelengths. For these two 
sources the spectral index and its error is calculated in the range 
between /190 cm and /13.6 cm. Point source no. 6 was resolved 
with VLA observati ons in B-configura tion into two components 
with 2" separation (Ca rilli et al.lll992| ). This source that has no 
known optical counterpart is likely a background AGN. 

4.2. Total flux density 

The total radio continuum flux density was integrated within el- 
lipses with the major axis aligned to the position angle (p. a. = 
52°) of NGC 253. The ellipses were chosen such that the visible 
extended emission in the disk and halo is covered. The parame- 
ters of the flux density integration are listed in Table [3] and the 
total flux densities obtained at different measurements are listed 
in Table|4] The flux densities of the point-like sources were sub- 
tracted from the total flux density, in order to obtain the flux 
density of the extended emission only. 

The average spectral index of the extended emission is a e = 
0.90 + 0.08 in the range between /190 cm and /13.6 cm. This 
spectrum is considerably steeper than that of young cosmic ray 
electrons (CREs) in supernova remnants (SNRs) with a « 0.5, 
although young s hell SNRs ca n have a steeper spectrum with 
a » 0.7 (see e.g. lGreenll2002l) . This is an indication for syn- 
chrotron losses of the CREs, which is be discussed in Sect. [53] 

Assuming a thermal fraction of 8 % at A20 cm dNiklas et alj 

we find a synchrotron non-thermal spectral index of a nt = 

1.0. 



4.3. Extended emission 

In this section we present the spectral index distribution of the 
extended emission, which we obtained from the total power 
maps. The center of all spectral index maps is dominated by the 
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Table 2. Unresolved sources which are denoted in Figs.[4]|7] The error of the flux densities was calculated assuming a 5 % calibration 
error. 



No. 


RA 


Dec 


5,i90 [mJy] 


Sm [mJy] 


5 ,16.2 


[mJy] 


S,\3.( 


[mJy] 


a 




1 


00 h 47 m 00 s .02 


-25°23'50".l 




4.1 ±0.2 


2.2 ± 


0.1 






0.5 ± 


0.2 


2 


00 h 47 m 12 s .01 


-25°17'43".9 




17.3 ± 1.0 


6.6 ± 


0.3 






0.8 ± 


0.2 


3 


00 h 47 m 33U2 


-25°17'17",3 


3960 ± 80 


2040 ± 100 


1270 


±60 


980 


±50 


0.44 


±0.02 


4 


00 h 47 m 44 s .91 


-25°13'38".4 




26 ±2 


7.9 ± 


0.4 






1.0 ± 


0.2 


5 


00 h 47 m 52 s .41 


-25°15'44",8 






1.5 ± 


0.2 










6 


00 h 47 m 58 s .90 


-25°18'31".9 


190 ± 10 


53 + 3 


19 ± 


1 


9.5: 


t0.5 


0.84 


±0.02 


7 


00 h 48 m 00 s .84 


-25° 16' 26". 3 




2.1 ±0.3 


1.10 


±0.06 






0.5 ± 


0.2 


8 


OO h 48 m 05 s .O 


-25°11'11".8 




5.3 ±0.3 


3.05 


±0.15 






0.4 ± 


0.2 



Table 3. The total flux densities were integrated within ellipses 
chosen to include the extended emission in the disk and halo. 
The major axis a, and the minor axis b for each wavelength A is 
presented. 



A [cm] 


a [arcsec] 


a [kpc] 


b [arcsec] 


b [kpc] 


90 


620 


12.1 


392 


7.6 


20 


620 


12.1 


392 


7.6 


6.2 


680 


13.3 


430 


8.4 


3.6 


648 


12.6 


400 


7.8 



Table 4. Integrated total power flux densities. Si is the total in- 
tegrated flux density and S e is the flux density of the extended 
emission only. The error of the flux densities was calculated as 
the quadratic sum of a 5 % calibration error and the baselevel 
error. 



A [cm] Instrument S, [ly] 5 C [ly] 

90 VLA 16.5 ±1.9 12.5 ±1.9 

20 VLA 6.3 ±1.1 4.2 ±1.1 

6.2 VLA + Effelsberg 2.71 ±0.14 1.44 ±0.07 

3.6 Effelsberg 1.66 ±0.08 0.68 ± 0.04 



nuclear point-like source (no. 3) with its flat spectral index of 
a = 0.44. In Fig. [8] we present the /190cm//l20cm spectral in- 
dex map. The A20 cm map was convolved with a Gaussian to 
match the resolution of the /190 cm map with 70". At this low 
resolution the inner disk and the spiral arms are not resolved. 
We find that the spectral index gradually steepens with increas- 
ing distance from the major axis, with the gradient strongest in 
the inner part of the disk. There is a slight trend that indicates 
that the southwestern part possesses a flatter spectral index than 
the northeastern part. At the southern tip of the radio spur the 
point-like source no. 6 is recognized by its flatter spectral index 
of a = 0.84 compared with the spectral index in the halo. 

The /I20cm//16.2cm spectral index map presented in Fig. [9] 
with its higher resolution of 30" reveals more details. The spiral 
arms are well resolved and show a flatter spectral index than 
the inter-arm regions. In addition the inner disk has a notable 
plateau with a flat spectral index of 6' extent along the major 
axis, which is not so clearly seen in the other maps. We find 
regions of steep spectral index (a > 1) at the base of the radio 
spur and at the northeastern extension. This finding agrees with 
the /I90cm//120cm spectral index map presented above, which 
has a steep spectral index at these locations, too. 

The A90 cm//13.6 cm spectral index map presented in Fig.fTOl 
has the lowest angular resolution of 84" among our spectral 
index maps. The morphology of this map is similar to the 
A20 cm//l6.2 cm map even though smoothed out in comparison. 
The spectral index is flat in the inner disk and in the spiral arms 
with a < 0.7. Limited by the sensitivity of the /13.6 cm map 
we find only four small extensions with extra-planar spectral in- 
dices, which are steeper than in the disk. The exception is the ex- 
tension northwest of the nucleus, which contains the point-like 
source no. 2 with a flat spectral index of a = 0.68. 



5. Cosmic ray transport 

While the CRs are created and accelerated in the disk, the ob- 
servation of extra-planar radio continuum emission requires a 
transport of CRs from the disk into the halo. This transport is 
investigated in this section. 

5. 1 . The scaleheight of the thick radio disk 

With the observed inclination angle of ; = 78.5° + 0.5° dPencd 
1980 ) any emission profile taken perpendicular to the disk is a 
superposition of disk and halo emission. We created a template 
of the vertical emission profile by projecting the emission pro- 
file along the major axis onto the minor axis. Moreover, the pro- 
file was convolved to the resolution of the observations in order 
to account for the beam smearing. Thus, the effective beam size 
was determined by fitting a Gaussian to the resulting profile. Any 
emission profile perpendicular to the disk that does not exceed 
the effective beam size cannot be attributed to extra-planar emis- 
sion but is just the observed inclined emission of the disk. Hence, 
we arrived at effective beam sizes of 1.2', 1.2', 2.1', and 2.5' for 
A90 cm, A20 cm, A6.2 cm, and A3. 6 cm, respectively. 

For the following analysis of the extended emission we sub- 
tracted all point-like sources by fitting Gaussians to them. The 
emission profiles were obtained using a strip integration perpen- 
dicular to the major axis with a strip width of 1.5' (1.7 kpc) and 
a pixel separation of 0.2 5' (0.29 kpc). With the routine discussed 
in iDumke et alj d 1995b profiles consisting of two exponential 
functions were fitted to the measured profiles taking the effec- 
tive beam size into account. This yields a scaleheight of the thin 
and the thick radio disks for the northern and the southern hemi- 
sphere, respectively. The synchrotron scaleheight h syn is calcu- 
lated as the mean of the scaleheight of the thick radio disk in 
the northern and the southern hemisphere. In case we could not 
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Fig. 8. Spectral index distribution between A90 cm and A20 cm 
with 70" resolution. Contours are at 0.6, 0.8, 1.0, 1.2, and 1.4. 
The total power radio continuum maps were clipped at a level of 
4x the rms noise prior to the combination, resulting in a spectral 
index error of less than +0.3. 
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Fig. 10. Spectral index distribution between A6.2 cm and A3 .6 cm 
with 84" resolution. Contours are at 0.2, 0.4, 0.8, 1.2, 1.6, and 
2.0. Both total power radio continuum maps were clipped at a 
level of 7x the rms noise prior to the combination, resulting in a 
spectral index error of less than ±0.5. 
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Fig. 9. Spectral index distribution between /t20cm and /16.2 cm 
with 30" resolution. Contours are at 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 
and 1 .8. Both total power radio continuum maps were clipped at 
a level of 5x the rms noise prior to the combination, resulting in 
a spectral index error of less than +0.4. 



obtain a two-component fit, we took the one-component expo- 
nential scaleheight as the thick radio disk. 

The scaleheights are smallest in the center and increase 
along the major axis further away from the center. This is 
shown in Fig.[TT| where also the variations of the total mag- 



netic field strength and the CRE lifetime are presented. The 
error of the scaleheights is estimated by varying the choice of 
data points. Moreover, the difference between the northern and 
the southern scaleheights contributes to the error. The weighted 
mean scaleheight at A6.2 of the thin radio disk is 0.33' ± 0.05' 
(0.38 + 0.06 kpc). The thick radio disk at A6.2 cm has a scale- 
height of 1.5' ± 0.1' (1.7 + 0.1 kpc). At /120 cm scaleheights of 
0.38'+0.05' (0.44+0.06 kpc) for the thin radio disk and 1.5'+0.1' 
(1.7 + 0.1 kpc) for the thick radio disk are found. Note that the 
/120 cm profiles drop beyond 3' height above the disk due to 
missing large-scale structure. However, this has only a marginal 
effect on the scaleheights. 

The /16.2 cm and /120 cm maps have a sufficient signal-to- 
noise (S/N) ratio, so that we were able to obtain two-component 
exponential distributions, whereas the /13.6cm map does not 
have a sufficient S/N ratio so that we used only a one-component 
exponential distributions. At /13.6cm we found a mean scale- 
height of 0.5' ± 0.1' (0.6 ± 0.1 kpc) and at A90cm we found 
2.2' + 0.2' (2.5 + 0.2 kpc). The scaleheights at A3.6 cm are domi- 
nated by the thin disk, because the low signal does not allow the 
thick disk to be traced to sufficiently large distances. In contrast 
at A9Q cm the emission is dominated by the thick radio disk so 
that one-component exponential fits are sufficient. 

The average scaleheight of NGC 253 agrees well with those 
reported for other edge-on galaxies at A6 cm (see also Sect. [6). 
That the scaleheight of the thick radio disk is larger than the 
effective beam proves the finding of extra-planar emission and 
hence the existence of a radio halo. The scaleheights of the thick 
disk at A6.2 cm, A20 cm, and A9Q cm are smallest for the central 
strips where the underlying SFR is highest. This important ob- 
servation is the result of synchrotron loss and is discussed in the 
next sections. 
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Fig. 11. Scaleheight h syn as a function of the offset x along the major axis (northeast is positive and southwest is negative). The 
magnetic field strength Bo (in units of fiG on the right axis) and the electron lifetime f e (in units of 10 6 yr on the right axis) are also 
shown. 



5.2. The total magnetic field strength 

From the non-thermal total power emission one can determine 
the total magnetic field strength assuming equipartition between 
the energy density of the CRs and the magnetic field. Here we 
used the revised equipartition formula bv lBeck & Krausej d2005l) 
assuming a pathlength of 6.5 kpc, which is the full width to the 
half power of the Gaussian fitted to the total power emission 
along the major axis. Moreover, we assumed a non-thermal spec- 
tral index of a nt = 1 .0 (see Sect. I4.21 i. a fractional polarization 
of 5.5 % (see Paper II), and a number density ratio Ko = 100 of 
protons and electrons. We stress that the last assumption might 
not be well fulfilled at large distances from the disk, because the 
CREs suffer high synchrotron losses whereas the CR protons are 
almost unaffected. Hence, the observed K factor increases with 
increasing distance from the CR sources which are located in 
the disk. The magnetic fiel d strength is then unde restimated by 
a small factor of (K/K ) l/4 dBeck & Krausdl2005l) . 

The structure of the regular magnetic field in NGC 253 is a 
superposition of a disk and a vertical component (see Paper II). 
The disk component is a spiral with a constant pitch angle and 
the inclination angle of the optical disk. The vertical component 
is orientated perpendicular to the disk. The above presented pa- 
rameters as well as the geometry of the magnetic field were as- 
sumed to be constant when computing the total magnetic field 
strength. To check the sensitivity of this assumption we varied 



the parameters in a certain range and observed its influence on 
the total magnetic field strength. We used an integration path 
length between 5 kpc and 20 kpc, a fractional polarization be- 
tween 5 % and 20 %, a non-thermal spectral index between 0.8 
and 1 .2, and assumed two structures for the magnetic field: disk- 
parallel or vertical. The former one describes the geometry of 
the magnetic field structure in the disk, where all possible mag- 
netic field directions along the line of sight are taken into ac- 
count. The latter one describes the geometry of the magnetic 
field in the halo, where a constant magnetic field direction along 
the line of sight is assumed. The variation of the total magnetic 
field strength in this range of conditions was found to be 1 .5 jjG, 
which we thus adopt as an estimate for the systematic error of 
the total magnetic field strength. 

The total magnetic field strength was scaled with the non- 
thermal s ynchrotron intensity 7 n using the standard equipartition 



tnermal synchrotron intensity l n i 
formula (Beck & Krause 2005): 



Beef 



l/(3+a nl ) 



(1) 



The distribution of the total magnetic field strength shown in 
Fig-EJwas calculated from the A6.2 cm map. It ranges between 
7 fiG and 18 fiG in the disk and is smaller in the halo. Note the 
subtracted nucleus, which is highlighted by a notable residual 
hole in the center of the map. The strip integration as described 
above applied to the map of the total magnetic field strength was 
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Fig. 12. Total power emission profiles in yuJy/beam at /16.2 cm 
perpendicular to the disk. The solid line shows the measured pro- 
file and the dashed line shows the fitted profile consisting of two 
exponential functions convolved with the effective beam. The 
offset of the profiles along the major axis is denoted. 



used in order to find the maximum of the total magnetic field 
strength Bo m eacn strip. This gives a lower limit on the lifetime 
of the CREs as described in the next section. 

IThompson et al.1 d2006l) argued that the magnetic field 
strength in starburst galaxies may be underestimated using the 
equipartition value. Our measured mean magnetic field strength 
of 12//G in the disk is accurately in equipartiti on with the 
mean gas surface-density 2 g = 0.84 X 10~ 2 gctrT 2 dPuche et al.l 
1991 s). so that the disk of NGC 253 fits onto the equiparti- 
tion line, although its starburst nucleus does not (see Fig. 1 in 
IThompson et al.l (120061) ). 



5.3. Loss processes of the cosmic ray electrons 

Several processes are known through which the CREs lose en- 
ergy, like synchrotron radiation, inverse Compton (IC) losses, 
non-thermal bremsstrahlung, ionization losses, and adiabatic 
losses. They all have different timescales which determine the 
lifetimes of the CREs. In order to calculate the timescales we 



adopt the equations bv lPohll ( Il990h . The synchrotron lifetime of 
CREs is 



f syn = 8.352 x 10 9 £(GeVr'fi(juGr 2 yr, 



(2) 



where E is the CRE energy and B is the total magnetic field 
strength. In the following we use the relation between the CRE 
energy E, the o bservation frequency v, and the total magnetic 
field strength B dRvbicki & Lightmanlll986l) : 



£(GeV) = (v/(16.1 MHz)) 1/2 fi(^G) 



-1/2 



(3) 



For the /16.2 cm observations with v = 4860 MHz and a typi- 
cal magnetic field strength of 15 fiG we find E = 4.5 GeV and 
hence f syn = 8.2 x 10 6 yr. Inverse Compton (IC) radiation losses 
have the same energy dependence as synchrotron losses and the 
corresponding lifetime is 



fic = 3.55 x 10 8 (£(GeV))- 1 (f/ ph /(10- 12 ergcm- 3 )r 1 yr, 



(4) 



where J/ P h is the tota l photon energy densi l 
frared observations of iRadovich et al.l (12001 1) 



y. From the in- 
we find a total 

photon energy density of 1.3 x 10~ 12 erg cirT 3 which leads to 
tic = 6.1 x 10 7 yr. This is almost an order of magnitude larger 
than the synchrotron lifetime. 

For the non-thermal CRE bremsstrahlung lifetime we use 



'br> 



3.96 x lOVcm"" 3 )" 1 yr, 



(5) 



where n is the ISM gas density. The ionization losses lead to a 
CRE lifetime of: 



? ion = 9.5 x 10 6 «(cm- 3 )- 1 £(GeV)yr. 



(6) 



We now use a mean gas colu mn density of £ e = 4 M Q pc~ 2 for 
the total gas as determined by iPuche et all d!991l) from H I ob- 
servations. We define a cylinder with the diameter of the disk 
D25 = 27.7' (31.7 kpc) and adopt a height of 3.4 kpc which is 
twice the scaleheight of the thick radio disk. The mean den- 
sity of p = 7.9 x 10~ 26 gcrrT 3 or n « 0.05 cm -3 leads to a 
bremsstrahlung lifetime of fbrems = 8 x 10 s yr and an ionization 
lifetime of t lon = 8.6 x 10 8 yr. 
The adiabatic timescale 



5I S) 



(7) 



depends only on the reciprocal velocity gradient. At this point we 
do not have any information about the velocity profiles. We esti- 
mate, however, the timescale of the adiabatic losses in Sect. 15.51 
as 3.6 x 10 7 yr, so that the synchrotron losses and the adiabatic 
losses are the most important loss processes of the CREs. 

The local synchrotron lifetime calculated with Eq.[2]may be 
an overestimate, because we observe at a certain location the su- 
perposition of locally produced CRs and CRs, which were trans- 
ported there from their origins. The latter ones have already lost 
a fraction of their energy and have shorter lifetimes than indi- 
cated by the local synchrotron losses. This effect is large in the 
halo, where almost all CRs were transported away from their ori- 
gins. In the disk, where we observe the highest fraction of locally 
produced CRs, this effect is small. Hence, we used the minimum 
of the local synchrotron lifetime f syn o calculated from the maxi- 
mum magnetic field strength Bo in the disk as the most accurate 
estimate. 

Because B decreases with z, f syn ,o is an underestimate. An 
improved estimate needs integration over the B profile, which is, 
however, not known well enough. As a starting point we assume 
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Fig. 13. Map of the total magnetic field strength with a resolution 
of 30". Contours are at 8, 12, 16, 20, and 24 x 1 /uG. 
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Fig. 14. Scaleheight h %yn of the thick radio disk as a function of 
the CRE lifetime ? e in the northeastern halo (between -6 kpc and 
kpc offset along the major axis). The linear fit is the theoretical 
expectation for a convective CR transport with a constant CR 
bulk speed. 



an exponential vertical profile. As shown in Appendix|A]we find 
in this case that r syn o is overestimated by no more than 10% if 
we restrict the data to z < h syn . 

Summarizing this section, the CRE energy loss is dominated 
by the synchrotron radiation and the adiabatic losses. Hence, we 
have to combine the two energy losses in order to estimate the 
lifetime f e of the CREs, using: 



1 1 1 

^syn,0 ^ad 



(8) a 



Since the adiabatic time scale is constant, the CRE lifetime is a 
function of the synchrotron lifetime. This can be seen in Fig.fTTl 
where we show the variation of Bo and f e along the major axis. 
The magnetic field strength has its maximum in the center of the 
disk where both the synchrotron lifetime and the CRE lifetime 
have a minimum. 



5.4. The cosmic ray bulk speed 

As described in Sect. [3] the total power emission resembles a 
dumbbell shape at all observed wavelengths; the contour lines 
show a strong thinning in the center close to the nucleus at the 
brightest part of the disk. The analysis of the scaleheight of 
the thick radio disk confirms this morphological analysis: it in- 
creases with galactocentric radius with the minimum at the strip 
centered on the nucleus. The opposite behavior is found for the 
magnetic field strength, which decreases from the center to the 
outer part of the disk. Hence, the scaleheight is small where the 
synchrotron lifetime f syn o and hence the CRE lifetime f e is small 
and vice versa. 

This relation motivates to calculate the average CR bulk 
speed, 



h e 

v = — , 

fe 
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Fig. 15. Same as Fig. [14] in the southwestern halo (between 
1.5 kpc and 6 kpc offset along the major axis). The linear fit 
(dashed line) is the theoretical expectation for a convective CR 
transport with a constant CR bulk speed, while the square root fit 
(solid line) is for a diffusive CR transport (see text for details). 



where h e is the CRE scaleheight. Inserting the non-thermal syn- 
chrotron intensity 7 n oc exp(-z//z syn ) into Eq. Q]and using the 
equipartition condition n e oc B 2 we find for the exponential pro- 
file of the CRE number density 



(9) n e °^exp|- — 

"syn 



2 

T+aZ 



(10) 
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where h syn is the synchrotron scaleheight of the thick radio disk. 
Thus, the CRE scaleheight is: 



h - 3 + a « h 

n e 2 "syn- 



(ID 



Since a nt = 1 the CRE scaleheight is twice the sclaleheight of 
the thick radio disk. 

In Fig. [14]the scaleheight of the thick radio disk at /190 cm, 
/120 cm, and /16.2 cm is shown for each strip as a function of 
the CRE lifetime. We find a (roughly) linear dependence be- 
tween /j syn and ? e . The CR bulk speed in the northeastern halo 
is 300 + 30kms _1 , where the error is dominated by the uncer- 
tainly in the magnetic field strength and thus in the electron 
lifetime. The linear least square fit has a reduced x 1 = 3.6. 
In the southwestern halo, we find a flatter dependence between 
the scaleheight and the CRE lifetime (Fig. [T5]) with a CR bulk 
speed of 250 ± 30kms _1 with a (x 1 = 2.7). In this part of the 
halo the scaleheight can be better fitted by square root function 

Or 2 = 0.8). 

In the northeastern halo, the CR bulk speed is constant be- 
tween A90 cm and A6.2 cm and does not vary with the CRE life- 
time. This can be explained by a convective CR transport, where 
the CRs and the magnetized ISM are transported together. In 
this case, the CR bulk speed is independent of the CR energy 
and hence in the northeastern halo the CR transport is probably 
dominated by convection. 

In the southwestern halo, the CR bulk speed decreases with 
increasing CRE lifetime and increasing wavelength. This can be 
explained by a diffusive CR transport where the CR bulk speed is 
v « k/(T c ■ h e ). Here, k is the diffusion coefficient and T c = 4/3 
is the adiabatic index of the relativistic CR gas. Replacing the 
electron scaleheight by h e — v ■ f e we obtain v * V/c/(r c • t e ). 
Since the diffusion coefficient is usually assumed to increase 
with the CR energy (k oc E° 5 ), the CR bulk speed decreases both 
with increasing wavelength and with increasing synchrotron life- 
time. This agrees with our observations and hence in the south- 
western halo diffusion probably plays a significant role. Using 
the square root fit of Fig. [15] we find a diffusion coefficient of 
k = (2.0 + 0.2) x 10 29 cm 2 s -1 which is in the range of the values 
quoted bv lBreitschwerdt et al.l d!993l) . 

The relation between the scaleheight and the CRE lifetime 
explains the observed pattern of the dumbbell total power emis- 
sion. In the center of the galaxy the magnetic field strength is 
highest and hence the synchrotron lifetime is smallest; it ham- 
pers the transport of CREs into the halo. This effect even over- 
comes the excess of emission in the center. The contour lines 
thus show a steep gradient in the center above and below the 
nucleus. 



5.5. Adiabatic losses 

Using the information about the CR transport we can now cal- 
culate the adiabatic loss timescale mentioned in Sect. 15.31 The 
CRs are accelerated to their bulk speed of 300kms _1 within 
the CRE scaleheight of 3.6 kpc. The adiabatic timescale is thus 
f a d = 3.6 x 10 7 yr using Eq. [7] assuming a constant accelera- 
tion. This timescale is similar to the synchrotron lifetime of the 
CREs. Hence, the adiabatic losses have to be considered in par- 
ticular for the low-energy CREs emitting at A90 cm. A transition 
of adiabatic loss dominated CR transport to radiative dominated 
CR transport should be visible in the spectrum as a steepening 
from a flat spectral index at low frequencies to a ste ep spectral 
index at high frequencies dPohl & Schlickeiserll990l) . The break 
in the spectrum with a steepening of Aa = 0.5 occurs at the 
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Fig. 16. Spectra of the diffuse emission in the northeastern halo, 
the central halo, and the southwestern halo. 



frequency where the adiabatic losses are equal to the radiation 
losses. 

Since the CR transport is different for the northeastern and 
southwestern halo it is worthwhile to investigate the spectra of 
the extended flux density in these regions separately. For this 
we integrated the flux density in three boxes with a width of 
5.2 kpc and a height of 5.7 kpc and an offset of -5.2 kpc, kpc, 
and 5.2 kpc along the major axis. The corresponding spectra 
in Fig. [16] steepen to higher frequencies in particular the spec- 
trum of the central region. This indicates high adiabatic losses 
at /190 cm. From linear fits we derive the spectral index as 
0.97 ± 0.05, 0.82 ± 0.05, and 0.89 + 0.05 in the northeastern, 
central, and southwestern halos, respectively. 

5.6. Cosmic ray spectral aging 

In the previous section we calculated the CR bulk speed from 
the scaleheight of the CRE and their lifetime. Since the radia- 
tion power of the synchrotron emission is proportional to E 2 , the 
spectral index steepens with time. This so-called CR spectral ag- 
ing can be used for an alternative approach to calculate the CR 
bulk speed as we show in this section. 

A single CRE loses its energy via synchrotron radiation with 
a loss power of 



dE _ o- T cBl i E \- 
dt 87r \mc 2 ) 



(12) 



where crj — 6.65 x 10~ 25 cm 2 is the Thomson cross section, Bo is 
the (constant) total magnetic field strength, and m — 51 1 keV/c 2 
is the electron rest mass. The time dependence of the CRE en- 
ergy is E{t) = Eo/(l + f/f syn ), where the synchrotron lifetime f syn 
is the time where a CRE has lost half of its initial energy Eq. 
The electron spectral index is defined by 



i Mi Ei 

r = -iog- log-, 



(13) 



where the number densities N\ and N2 are calculated for two dif- 
ferent CR energies E\ and £V The electrons are injected at t — 
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Fig. 17. Profiles of the spectral index between A90 cm and 
A6.2 cm perpendicular to the disk as a function of distance from 
the major axis. The offset of the strips along the major axis is 
denoted. 



and are assumed not to be reaccelerated in the halo. If we neglect 
diffusion and assume a constant wind speed, the electrons can be 
considered to reside in a box while they are losing their energy 
due to synchrotron r adiation. In this case the ir number density 
can be expressed by dLisenfeld & Volkll2000l) 



N(E) = K E- p 1 



(14) 



where p is the injection spectral index which is assumed to be 
about p ss 2. Our model describes the geometry of a "disk 
wi nd" with a vertical wind d irection. Our analysis is similar 
to ISeaquist & Odegardl d 1 99 lb who made a spherical symmet- 
ric model for the wind in the central part of M82. For the time 
derivative of the electron spectral index we find for t — > 



Ay /l 1\ / E 2 



(15) 



where t\ and t% denote synchrotron lifetimes. The radio spectral 
index can be computed using a = (y — l)/2. With Eqs.[2]and[3] 



we can express the time derivative of the radio spectral index by 

Aa _ BqQiG) 3 ' 2 

~At ~ ln(vi/v 2 )- 8.352 x 10 9 yr 



16.1MHz 



16.1MHz 



(16) 



which only depends on the total magnetic field strength and the 
observation frequencies v\ and v 2 . Thus, we get an expression 
for the average CR bulk speed v = Az/Af 



Aa / Aa 
~At Az' 



(17) 



where Aa/Az is the slope of the spectral index profiles. 

We applied a strip integration in order to calculate profiles of 
the spectral index between A9Q cm and A6.2 cm. The strip width 
was 180" and the pixel separation was 15". The spectral index 
profiles are presented in Fig.[T7]where the error bars were com- 
puted from the noise in the strips, which is smaller than the noise 
in the map (since the width of the strips is larger than the map 
resolution of 70"). The spectral index distribution shows an al- 
most linear steepening in some (but not all) of the strips. This 
can be explained by a constant time derivative Aa/At and hence 
a constant CR bulk speed v. 

The central strips at Okpc offset from the nucleus show a 
very flat spectral index for z < 1 kpc. This can be explained by a 
thermal flattening in the nuclear region at A90 cm dCarillilll996l ). 
We restricted the analysis of the spectral index slope to heights 
larger than 1 .4 kpc where we do not expect any contribution of 
the thermal gas in the disk. The resulting CR bulk speeds are 
shown in Table [5] where the subscripts "N" and "S" denote the 
northern and and southern hemisphere, respectively. The mean 
CR bulk speed derived from the spectral index slope is 170 ± 
70 km s , to be compared with the CR bulk speed obtained from 
the scaleheights of v,\(,.2 = (300 + 30)kms _1 . Since the time 
derivative is taken at t — this is a lower limit for the CR bulk 
speed. 



Table 5. Average CR bulk speeds obtained from the radio spec- 
tral index aging. 



X 

[kpc] 


Aal 
Az In 
[kpc- 1 ] 


Aa\ 
Az Is 
[kpc- 1 ] 


Bo 
IpG] 


Vn 

[km s- 1 ] 


vs 

[km s- 1 ] 


-6.9 


-0.22 


-0.12 


11 


107 


196 


-3.4 


-0.26 


-0.31 


14 


130 


109 





-0.30 


-0.26 


17 


150 


174 


3.4 


-0.16 


-0.12 


15 


234 


312 


6.9 


-0.10 


-0.20 


10 


204 


102 



6. Discussion 

The CR distribution to constrain the transport of CRs from the 
disk into the halo. The local CRE lifetime is dominated by syn- 
chrotron losses and depends on the magnetic field strength in the 
disk which is highest in the central part. The dependence of the 
scaleheight on the CRE lifetime, as evident from the dumbbell 
shaped halo, requires a vertical CR transport. 

The CR bulk speed, defined by the ratio of the CRE scale- 
height to their lifetime, is constant in the northeastern halo with 
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300 ± 30kms _I and remarkably similar when applied to our 
observations at three different wavelengths. Thus, the local CR 
bulk speed is independent of the CR energy, the star-formation 
rate (SFR), and the magnetic field strength. Our second method 
to determine the CR bulk speed is to use the spectral aging of 
the CREs. The steepening of the observed radio spectral in- 
dex with increasing distance from the galactic midplane results 
in a lower limit for the CR bulk speed with 170 + 70kms~ I . 
IZirakashvili & Volkld2006l) made an analytical model for the CR 
transport in the nuclear outflow. They found by comparison with 
radio continuum data a velocity of 300kms~' near the disk and 
900kms _1 as terminal velocity. 

iDahlem et all d 19951) and iDahlem et ail d2006l) studied the 
CR distribution in edge-on galaxies. They found a tight and lin- 
ear correlation between the radial extent of radio halos and the 
level of star-formation activity in the underlying disk. However, 
radio halos can be better studied by their vertical intensity dis- 
tributions. Whereas the amplitude of the vertical profile depends 
on the overall radio intensity of the disk, related to the SFR, its 
scaleheight is determined by the outflow velocity and the CRE 
lifetime. The observed vertical extent also depends on the signal- 
to-noise ratio of the measurements. Only the scaleheight is suited 
to characterize the formation of radio halos. 

Global radio scaleheights observed at A6 cm in three edge- 
on galaxies with very different global SFRs were found to be 
similar, about 1 .8 kp c for the thick radio disk (iDumke & Krausd 
[19981: iKrau se 2004). This agrees well with our result of 1 .7 kpc 
of NGC 253 which has the highest SFR of this sample. Our Eq.|9] 
indicates that a constant CRE scaleheight h e in galaxies with dif- 
ferent SFRs and hence different CRE lifetimes f e can only be 
achieved if the global CR bulk speed v increases with global 
SFR. As all four galaxies show a dumbbell shaped radio halo, the 
local scaleheight depends on the local magnetic field strength. 

Since we postulate an outflow, we analyse the pressure con- 
tributions of the different ISM phases. These are the thermal gas, 
the kinetic energy, and the magnetic field. The thermal pressure 
of the halo gas is Ptherm = InJi^T — 7 x 10~ 12 erg cm -3 , where 
we used an electron densi ty of n e — 8 x 10 ~ 3 c m" 3 and a tem- 
perature of T = 3 x 10 6 K (IStrickland et al.l2002l) . The ram pres- 
sure of the outflow is P rdm = l/2« e m p v 2 = 6 x 10~ 12 erg cm -3 , 
where ra p is the proton mass and v = 300 km s the wind veloc- 
ity. Assuming a typical magnetic field strength of B — 15 yuG, we 
find a magnetic pressure of P mag = B 2 /(8n) = 9xl0~ 12 ergcm" 3 . 

If the gas in the disk were in a hydrostatic equilibrium, the 
gas pressure would be P erav ~ 7rG£gX t - We obtain P grav = 1.3 X 
10" ergcirT 3 , where we inserted E g = 4M pc~ 2 for the gas 
surface-density dPuche et all 19911) . and£m = 7.4 4xl0~ 2 g cm -2 
for the total mass surface-density dPence 11981b . We note, that 
the thermal pressure, the ram pressure, and the magnetic pres- 
sure are of the same order of magnitude. Since the sum of the 
thermal pressure, the ram pressure, and the magnetic pressure 
is 2.2 x 10~" erg cm" 3 , this may indicate that the disk is not in 
hydrostatic equilibrium and allows outflow. 

If the outflow traced by the CR transport is a galactic wind, 
it should finally leave the gravitational potential of the galaxy. In 
this case the C R bulk speed must e xceed the escape velocity of 
v eS c — 2^ 2 -Qr dVeilleux et alj2005b . where fir is the rotation ve- 
locity of the disk. Of course, the outflow may accelerate further 
away from the disk, below our sensitivity limit, and thus a galac- 
tic wind is possible even if the observed CR bulk speed is too 
low. On the other hand, a CR bulk speed above the escape veloc- 
ity supports the galactic wind scenario, because it will certainly 
escape from the gravitational potential, unless the outflow decel- 



erates in the halo. Us ing the rotation speed of 200 km s" 1 mea- 
sured bv lPencel(ll981l) . we find an escape velocity of 280 km s _1 , 
which agrees within the errors with our measured CR bulk speed. 
Hence, a galactic wind is in NGC 253 likely and the CR transport 
shows the presence of a "disk wind" in the northeastern halo. 

The disk wind can efficiently transport gas into the halo. The 
transport should be more efficient in the convective northeastern 
halo. Indeed, there is significantly more neutr al hydrogen gas 
in the northeastern halo, as HI observations by iBoomsma et al] 
(120051) show. A similar distribution is found for the ionized hy- 
drogen, wher e prominent Hg plum es are found only in the north- 
eastern halo dHoopes et alj|1996l) . This picture is corroborated 
by huge lobes of soft X-ray emission in t he northeastern halo , 
which are absent in the southwestern halo dPietsch et al.ll2006T) . 
The discrepancy between the northeastern and the southwestern 
halo may be explained by the level of star-formation activity in 
the disk. Both the Ha emission and radio continuum emission 
is stronger in the northeastern part of the disk. This hints to a 
higher SFR that can more easily drive the disk wind against the 
gravitational pull. 

The disk wind may finally escape from the gravitational po- 
tential to the intergalactic medium (IGM). From the observations 
we cannot make a definite conclusion because we do not observe 
the transition between the radio halo and the IGM further away 
from the detected radio halo. With our observations we are not 
able to conclude whether the disk wind accelerates or deceler- 
ates within distance from the disk. With low frequency observa- 
tions the scaleheight of the synchrotron emission is expected to 
be larger and hence the radio halo could be observed to larger 
vertical heights. Such data could allow us to determine whether 
the disk wind accelerates as expected for a CR driven wind. 
Although LOFAR observations would be ideal for this purpose, 
the low declination of NGC 253 (-25°) requires to use the EVLA 
or the future SKA at frequencies as low as possible. 

7. Conclusions 

Our radio continuum study of NGC 253 shows an extended ra- 
dio halo at the wavelengths from /190 cm to /T3. 6 cm. Because 
studies of NGC 253 need a high dynamic range, it is a difficult 
target for radio continuum observations. We treated the high dy- 
namic range with a specially tailored data reduction technique. 
Our maps presented here are thus essentially noise limited. The 
combination of our new VLA mosaic map with the single-dish 
Effelsberg data in order to fill in the missing zero-spacing flux 
density allowed us to study the extra-planar continuum emission 
with high resolution and sensitivity. They are used to put de- 
tailed constraints on the CR transport, which lead to our main 
conclusions: 

1 . The scaleheight of the thick radio disk in the northeastern 
halo is a linear function of the CRE lifetime in the under- 
lying disk at /190 cm, /120 cm, and /16.2 cm. This can be un- 
derstood if the CR transport proceeds preferentially in the 
vertical direction. 

2. The steepening of the spectral index, which is a linear func- 
tion of z, can be explained by strong synchrotron losses 
of the CREs. The contribution of ionization losses and 
bremsstrahlung losses can be neglected. The adiabatic losses 
are of the same order of magnitude as the synchrotron losses 
at A20 cm but dominate at A9Q cm. 

3. We computed the CR bulk speed as the ratio of the total 
power scaleheight and the CRE lifetime. The CR bulk speed 
in the northeastern halo of 300+30 km s is similar to the es- 
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cape velocity of 280 km s , indicating the existence of a disk 
wind which escapes from the gravitational potential. The lo- 
cal CR bulk speed is independent of the CR energy, the SFR, 
and the magnetic field strength. 

4. The CR transport is mainly convective in the northeastern 
halo while in the southwestern halo the diffusion of CRs is 
significant. We measure an average diffusion coefficient of 
2.0 + 0.3 x 10 29 cm 2 s- 1 

5. The disk wind is the driving source for the luminous gas in 
the halo. In the convective northeastern halo the gas is trans- 
ported more efficiently from the disk into the halo than in 
the diffusive southwestern halo. This can explain the differ- 
ent amounts of extra-planar gas in these two halo parts as 
observed in H I, Ha, and soft X-rays. 

The CR transport is always the superposition of convection and 
diffusion. The data presented here provide a basis for more de- 
tailed calculations using the combined diffusion and convection 
equation. These calculations should include vertical profiles of 
the magnetic field strength, the CR number density, and the ve- 
locity of the galactic wind. 
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Appendix A: The synchrotron lifetime within an 
exponential magnetic field profile 



In order to calculate the CR bulk speed in Sect. l5.4l we used the 
synchrotron lifetime of the electrons. Since we do not know the 
vertical distribution of the magnetic field, we assumed a constant 
magnetic field strength. In reality, however, we expect the mag- 
netic field strength to have a vertical, decreasing profile. In this 
section we investigate the difference between the synchrotron 
lifetime for a constant magnetic field and for an exponential ver- 
tical profile. 

The exponential profile of the magnetic field reads 



B(z) = B Q exp(-z/h B ), 



(A.l) 



where Bq is the maximum magnetic field strength and h B the 
scaleheight of the magnetic field. Replacing t = z/v, where we 
assume a constant CR bulk speed, the differential of the energy 
is: 



dE = - 



cr T cB(zf 



E(z)Y dz 
mc 2 J v 



(A.2) 



We integrate this differential equation by separating the variables 
with: 



dE cr T B(z) 2 



dz. 



E 2 87to 2 c 3 v 
Now we integrate from z = to 
-E JI7 rz o-^Bl 



}§-f 

Je 11 Jz=0 



87rm 2 c 3 v 



exp(-2z//z B )<fe, 



(A.3) 



(A.4) 
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where we inserted the exponential profile of B. Finally, we find 
for the vertical CRE energy profile 



E(z) 



EqIibo-jB}, 

5-^ ■ (1 - exp(-2z/h B )) 



(A.5) 



1 



167rm c v 



If we use the definition for the synchrotron lifetime for a constant 
field strength Bq 



'syn,0 — 



87rm 2 c 3 



we find 



Eb(z) 



1 + 



2f S yn,0V 



(1 - exp(-2z//z B )) 



(A.6) 



(A.7) 



For the case of a constant magnetic field strength we have as 
energy profile: 



Ebo(z) 



Eo 



1 + 



(A.8) 



fsyn.OV 



Clearly, the larger the scaleheight of the magnetic field, the 
smaller is the difference between the two CRE energy profiles. 
We now use Iib = 4 ■ h syB as a lower limit due to the equipartition 
condition (see Eq.Q]), where we set a at = 1. In Fig. [AT] we plot 
the ratio of the energies Ebo/Eb as a function of z in units of 
Kb- Any estimate of the error made using a constant magnetic 
field instead of a exponential magnetic field must be averaged 
using the weight of the points. If we restrict z < h syn and take 
into account that the inner points have larger weights (due to the 
higher synchrotron emission), the error estimate is around 10 %. 
Hence, in our case the assumption of a constant magnetic field 
strength for the calculation of the CRE energy and hence of the 
synchrotron lifetime in Sect. 15 .3 l is acceptable. 
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Fig.A.l. Ratio of the two energies Ebo/Eb as a function of z. 
For Eb the magnetic field has an exponential distribution with a 
scaleheight lis while for £bo it is constant. 



